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Abstract 

Pulsars are among the most mysterious astrophysical objects in the Universe and 
are believed to be rotating neutron stars formed in supernova explosions. They are 
unique testing grounds of dense matter theories and gravitational physics and also 
provide links among nuclear physics, particle physics and General Relativity. Neutron 
stars may exhibit some of the most extreme and exotic characteristics that could not 
be found elsewhere in the Universe. Their properties are largely determined by the 
equation of state (EOS) of neutron-rich matter, which is the chief ingredient in calcu- 
lating neutron star structure and properties of related phenomena, such as gravitational 
wave emission from deformed pulsars. Presently, the EOS of neutron-rich matter is 
still very uncertain mainly due to the poorly known density dependence of the nuclear 
symmetry energy especially at supra-saturation densities. Nevertheless, significant 
progress has been made recently in constraining the density dependence of the nuclear 
symmetry energy mostly at sub-saturation densities using terrestrial nuclear reactions. 
While there are still some uncertainties especially at supra-saturation densities, these 
constraints could provide useful information on the limits of the global properties of 
pulsars and the gravitational waves to be expected from them. Here we review our 
recent work on constraining properties of pulsars and gravitational radiation with data 
from terrestrial nuclear laboratories. 
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1 Introduction 



Pulsars exhibit a large array of extreme characteristics. They are generally accepted to 
be rotating neutron stars - the smallest and densest stars known to exist. Matter in their 
cores is compressed to huge densities ranging from the density of normal nuclear matter, 
po « 0.16/m^ 1 , to an order of magnitude higher [1]. The number of baryons forming a 
neutron star is in the order of A ps 10 57 . Understanding properties of matter under such 
extreme conditions of density (and pressure) is still far from complete and represents one 
of the most important but also challenging problems in modern physics. 

Because of their strong gravitational binding neutron stars can rotate very fast [2]. The 
first millisecond pulsar PSR1937+214, spinning at v = 6A2Hz 0, was discovered in 
1982, and during the next decade or so almost every year a new one was reported. In 
the recent years the situation changed considerably with the discovery of an anomalously 
large population of millisecond pulsars in globular clusters where the density of stars 
is roughly 1000 times that in the field of the galaxy and which are therefore very favorable 
sites for formation of rapidly rotating neutron stars which have been spun up by the means 
of mass accretion from a binary companion. Presently, the number of the observed pulsars 
is close to 2000, and the detection rate is rather high. 

In 2006 Hessels et al. @ reported the discovery of a very rapid pulsar J1748-2446ad, 
rotating at v = IVoliz and thus breaking the previous record (of 6A2Hz). Rapid rotation 
could affect significantly the neutron star structure and properties, especially if the rota- 
tional frequency is near the Kepler frequency for the star. (The Kepler, or mass-shedding 
frequency, is the highest possible frequency an object bound by gravity can have before it 
starts to loose mass at the equator.) Pulsars with masses above 1M enter the rapid-rotation 
regime, if their rotational frequencies are higher than 1000Hz [2]. Although, presently no 
pulsar is confirmed to rotate above the 1000Hz limit, theory does not exclude the possi- 
bility for existence of such rapidly rotating neutron stars and therefore it is important to 
predict their properties. (Here we should mention the recent discovery of X-ray burst os- 
cillations from the X-ray transient XTE J1739285 [6], which could suggest that it contains 
an extremely rapidly rotating neutron star spinning at 1122 Hz. However, this observation 
has not been confirmed yet.) Neutron stars are natural astrophysical laboratories of dense 
matter [4]. Rotating neutron stars appear as much better probes for the structure of dense 
baryonic matter than static ones, primarily because the particle compositions in rotating 
neutron stars are not frozen in, as it is the case for non-rotating neutron stars, but are vary- 
ing with time [7]. The associated density changes could be as large as 60% [4] in neutron 
stars in binary stellar systems (e.g., low-mass X-ray binaries), which are being spun up to 
higher rotational frequencies, or isolated rotating neutron stars (e.g. isolated millisecond 
pulsars) which are spinning down to lower frequencies because of the gravitational radia- 
tion, electromagnetic dipole radiation, and a wind of electron-positron pairs Q . 

Since neutron stars are objects of extremely condensed matter, the geometry of space- 
time is considerably altered from that of a flat space. Therefore, the construction of realistic 
models of neutron stars has to be done in the framework of General Relativity [ 1 , 4 ] . De- 
tailed knowledge of the equation of state (EOS) of stellar matter over a very wide range of 



4 



Plamen G. Krastev and Bao-An Li 



densities is required for solving the neutron-star structure equations. At present time the be- 
havior of matter under extreme densities such as those found in the interiors of neutron stars 
is still highly uncertain and relies upon, often, rather controversial theoretical predictions. 
On the other hand, fortunately, heavy-ion reactions provide unique means to constrain the 
EOS of dense nuclear matter in terrestrial laboratories [8]. One of the major uncertainties 
of the EOS is the density dependence of the nuclear symmetry energy |j9l [TOl [TTJ , E sym {p), 
which is the difference between the nucleon specific energies in pure neutron matter and 
symmetric nuclear matter. Due to its importance for the neutron star structure and many 
other ramifications in astrophysics and cosmology, determining the density dependence 
of the nuclear symmetry energy has been a major thrust of research for the intermediate 
energy heavy-ion community. Although extracting information about the density depen- 
dence of the nuclear symmetry energy from heavy-ion reactions is not an easy task due 
to the complicated roles of the isospin degree of freedom in the reaction dynamics, sev- 
eral promising probes of the symmetry energy have been suggested and tested in recent 
yearsIEl El El El El (see also Refs. 00 El COD for reviews). In fact, some signifi- 
cant progress has been made in determining the E sym (p) at subsaturation densities using: 
(1) isospin diffusion |[T9l and isoscaling [20] in heavy-ion reactions at intermediate ener- 
gies |2l]|22l|23ll2S|25l, (2) sizes of neutron skins in heavy nuclei flTOl E2l EH , (3) the 
Pygmy dipole resonance ||29ll and (4) masses of nuclei IT3011 . At supranormal densities, a 
number of potential probes of the symmetry energy have been proposed ll"i~6l [3~T1 although 
the available data has been very limited so far. Interestingly, circumstantial evidence for 
a super-soft E sym {p) at supra-saturation densities was reported recently [32] based on the 
analysis of the FOPI/GSI data on pion production |[33l . It is also very encouraging to no- 
tice that several dedicated experiments are being planned at several laboratories to study in 
more detail the symmetry energy at supra-saturation densities using high energy heavy-ion 
reactions induced by both stable and radioactive beams. 

While global properties of spherically symmetric static neutron stars have been studied 
extensively and comprehensive literature exists, e.g. Refs. ll9l fTOl ITT1 [34l l35l [36l l37l [38l . 
properties of (rapidly) rotating neutron stars have been investigated to lesser extent. Models 
of (rapidly) rotating neutron stars have been constructed by several research groups with 
various degree of approximation S|351SagDEl@31ElElEIEl^l^llSD (see 
Ref. lIBTl for a review). (Rapidly) rotating, elliptically deformed neutron stars could be 
one of the major candidates for emitting gravitational waves (GWs) - tiny ripples in space- 
time predicted by the theory of General Relativity 11521 . Although GWs have not been 
detected directly yet, indirect evidence do exist ll53l . GWs are characterized by a small 
dimensionless strain amplitude, ho, which, in addition to the pulsar's distance to detector, 
depends on the neutron star structure determined by the underlying EOS of stellar matter. 
Thus the EOS is instrumental for predicting the strength of gravitational radiation emitted 
from neutron stars. 

The partially constrained EOS by the available terrestrial nuclear laboratory data has 
been used in studying several properties of neutron stars, including the mass-radius corre- 
lation [54], the surface temperature of neutron stars in connection with the changing rate 
of the gravitational constant G [55], the core-crust transition density [56], the frequency 
and damping time of the w-mode of gravitational waves ll57l and the gravitational binding 
energy of neutron stars ll58l . In this chapter we concentrate on reviewing several selected 
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properties of fast pulsars and the gravitational waves expected from them mostly based on 
our recent work reported in Refs. l59ll60ll6lTl . 

2 The equation of state of neutron-rich nuclear matter partially 
constrained by recent data from terrestrial heavy-ion reac- 
tions 

In this section, we first outline the theoretical tools one uses to extract information about 
the EOS of neutron-rich nuclear matter from heavy-ion collisions. We put the special em- 
phasis on exploring the density-dependence of the symmetry energy as the study on the 
EOS of symmetric nuclear matter with heavy-ion reactions is better known to the astro- 
physical community and it has been extensively reviewed, see e.g., Refs. JU [lOl |62l for 
recent reviews. We will then summarize the latest constraints on the density dependence 
of the symmetry energy at sub-saturation densities extracted from heavy-ion reactions at 
intermediate energies. Finally, we address the question of what kind of isospin-asymmetry, 
especially for dense matter, can be reached in heavy-ion reactions. 

Heavy-ion reactions provide unique means to create dense nuclear matter in terrestrial 
laboratories similar to those found in the core of neutron stars. While unlike the matter 
in neutron stars the dense matter created in heavy-ion reactions is hot, accurate informa- 
tion about the EOS of cold matter can be extracted reliably by modeling carefully the 
kinetic part of the EOS during heavy-ion reactions. Depending on the beam energy, im- 
pact parameter and the reaction system, various hadrons and/or partons may be created 
during the reaction. To extract information about the EOS of dense matter from heavy- 
ion reactions requires careful modeling of the reaction dynamics and selection of sensitive 
observables. Among the available tools, computer simulations based on the Boltzmann- 
Uehling-Uhlenbeck (BUU) transport theory have been very useful, see, e.g., Refs. ||8l |63T| 
for reviews. The evolution of the phase space distribution function /j(r, p, t) of nucleon i is 
governed by both the mean field potential U and the collision integral I collision via the BUU 
equation 

dfi - 

-^T+VpU- V r fi - V r U ■ V p fi = I co lli sl on. (1) 

Normally, effects of the collision integral I co iiision y i a both elastic and inelastic channels 
including particle productions, such as pions, are modeled via Monte Carlo sampling using 
either free-space experimental data or calculated in-medium cross sections for the elemen- 
tary hadron-hadron scattering ll63l . The collision integral is critical for modeling the kinetic 
part of the EOS. Information about the EOS of cold nuclear matter is obtained from the un- 
derlying mean-field potential U which is an input to the transport model. By comparing 
experimental data on some carefully selected observables with transport model predictions 
using different mean-field potentials corresponding to various EOSs, one can then con- 
strain the corresponding EOS. The specific constrains on the density dependence of the 
nuclear symmetry energy that we are using in this work were obtained by analyzing the 
isospin diffusion data |[T9ll within the IBUU04 version of an isospin and momentum depen- 
dent transport model ll64l . In this model, an isospin and momentum-dependent interaction 
(MDI) ll65l is used. With this interaction, the potential energy density V(p, T, 5) at total 
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density p, temperature T and isospin asymmetry 5 is 

yp,r,5 = + ^(p n + p p )+ (i-xS) 

Po 2p p a + 1 pi 

In the mean field approximation, Eq. © leads to the following single particle potential for 
a nucleon with momentum p and isospin r 

U T (p,T,6,p,x) = A u (x)^ + Mx)^ + b(P-Y (l-x5 2 ) 

Po Po \PoJ 

a + l pi t ±^_ r po J l + (p-py/k l 

where r = 1/2 (—1/2) for neutrons (protons), x, A u {x), A^{x), B, C T:T ,C Tj - T , a, and 
A are all parameters given in Ref. [65]. The last two terms in Eq. © contain the mo- 
mentum dependence of the single-particle potential, including that of the symmetry po- 
tential if one allows for different interaction strength parameters C Tj _ T and C TjT for a 
nucleon of isospin r interacting, respectively, with unlike and like nucleons in the back- 
ground fields. It is worth mentioning that the nucleon isoscalar potential estimated from 
Ui SO scaiar ~ (U n + U p )/2 agrees with the prediction of variational many -body calculations 
for symmetric nuclear matter ll66l in a broad density and momentum range ll64l . Moreover, 
the EOS of symmetric nuclear matter for this interaction is consistent with that extracted 
from the available data on collective flow and particle production in relativistic heavy-ion 
collisions up to five times the normal nuclear matter [8]. On the other hand, the corre- 
sponding isovector (symmetry) potential can be estimated from U sym (U n — U p ) /25. At 
normal nuclear matter density, the MDI symmetry potential agrees very well with the Lane 
potential extracted from nucleon-nucleus and (n,p) charge exchange reactions available for 
nucleon kinetic energies up to about 100 MeV [64]. At abnormal densities and higher 
nucleon energies, however, there is no experimental constrain on the symmetry potential 
available at present. 

The different x values in the MDI interaction are introduced to vary the density de- 
pendence of the nuclear symmetry energy while keeping other properties of the nuclear 
equation of state fixed. Specifically, choosing the incompressibility Kq of cold symmet- 
ric nuclear matter at saturation density p$ to be 211 MeV leads to the dependence of the 
parameters A u and A[ on the x parameter according to 

, / x 25 2B 

A u (x) = -95.98 - x -, Ai(x) = -120.57 + x (4) 

<7 + l a + 1 

with B = 106.35 MeV. 

With the potential contribution from Eq. © and the well-known contribution from nu- 
cleon kinetic energies in the Fermi gas model, the EOS and the symmetry energy at zero 
temperature can be easily obtained. As shown in Fig. CO, varying the parameter x leads to 
a broad range of the density dependence of the nuclear symmetry energy, similar to those 
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Figure 1 : The density dependence of the nuclear symmetry energy for different values of 
the parameter x in the MDI interaction. Taken from f72l . 

predicted by various microscopic and/or phenomenological many-body theories. As previ- 
ously demonstrated in [24] and [54], only equations of state with x between -1 and have 
symmetry energies in the sub-saturation density region consistent with the isospin diffusion 
data and the available measurements of the skin thickness of 208 Pb using hadronic probes. 
More recently however, it was shown that at supra-saturation densities the symmetry energy 
could be much softer than the extrapolation of the symmetry energy with x = ll32l . The 
possibility for a super-soft symmetry energy at supra-saturation densities could have impor- 
tant consequences for the neutron star stability and the nature of gravity at sub-millimeter 
distances ll67l . Nevertheless, as the first and a conservative step in our studies of the prop- 
erties of (rotating) neutron stars ||59l 1601 and the gravitational waves expected from them 
EQ, we have used x = —I and x = to compute the range of possible (rotating) neutron 
star configurations. The E sym {p) with x = is also consistent with the RMF prediction 
using the FSUGold interaction [68]. In this review, we thus consider only the two limiting 
cases with x = and x = — 1 as boundaries of the symmetry energy consistent with the 
available terrestrial nuclear laboratory data at sub-saturation densities. 

To facilitate comparisons with other models in the literature, it is useful to parameterize 
the E sym (p) from the MDI interaction and list its characteristics. Within phenomenological 
models it is customary to separate the symmetry energy into the kinetic and potential parts, 

Table 1: The parameters F (MeV), G , K sym (MeV), L (MeV), and K asy (MeV) for different 
values of x. Taken from [22 1. 



X 


F 


G 


lv sym 


L 


K 


1 


107.232 


1.246 


-270.4 


16.4 


-368.8 





129.981 


1.059 


-88.6 


62.1 


-461.2 


-1 


3.673 


1.569 


94.1 


107.4 


-550.3 


-2 


-38.395 


1.416 


276.3 


153.0 


-641.7 
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Figure 2: Central baryon density (upper panel) and isospin asymmetry (lower panel) of high 
density region in the reaction of 132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon and 
an impact parameter of 1 fm. Taken from Ref . f72l . 

see, e.g. ||69ll , 

E sym {p) = (2 2 / 3 - lf-E% (p/ P0 ) 2/3 + E%tn(p). (5) 

With the MDI interaction, the potential part of the nuclear symmetry energy can be well 
parameterized by 

EfyL(p) = F(x)p/ P0 + (18.6 - F{x)){p/p,) G ^\ (6) 

with F{x) and G(x) given in Table[T]for x = 1, 0, —1 and —2. The MDI parameteriza- 
tions for the Efy^p) are similar to but significantly different from those used in Ref. ||69l . 
Also shown in Table Q] are other characteristics of the symmetry energy, including its slope 
parameter L and curvature parameter K sym at po, as well as the isospin-dependent part 
iv" asy of the isobaric incompressibility of asymmetric nuclear matter (22). The symmetry 
energy in the subsaturation density region with x=0 and - 1 can be roughly approximately 
by E sym (p) « 31.6(p/ / o ) - 69 and E sym (p) « 31.6(p/ / o ) 1 ' 05 , respectively. 

What is the maximum isospin asymmetry reached, especially in the supra-normal den- 
sity regions, in typical heavy-ion reactions? How does it depend on the symmetry energy? 
Do both the density and isospin asymmetry reached have to be high simultaneously in or- 
der to probe the symmetry energy at supra-normal densities with heavy-ion reactions? The 
answers to these questions are important for us to better understand the advantages and lim- 
itations of using heavy-ion reactions to probe the EOS of neutron-rich nuclear matter and 
properly evaluate their impacts on astrophysics. 

To answer these questions we first show in Fig.[2]the central baryon density (upper win- 
dow) and the average (n/p) p > po ratio (lower window) of all regions with baryon densities 
higher than po in the reaction of 132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon and 
an impact parameter of 1 fm. It is seen that the maximum baryon density is about 2 times 
the normal nuclear matter density. Moreover, the compression is rather insensitive to the 
symmetry energy because the latter is relatively small compared to the EOS of symmetric 
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nuclear matter around this density. The high density phase lasts for about 15 fm/c from 5 to 
20 fm/c for this reaction. It is interesting to see in the lower window that the isospin asym- 
metry of the high density region is quite sensitive to the density dependence of the symmetry 
energy used in the calculation. The soft (e.g., x = 1) symmetry energy leads to a signifi- 
cantly higher value of {n/p) p > Po than the stiff one (e.g., x = —2). This is consistent with 
the well-known isospin fractionation phenomenon in asymmetric nuclear matter iTTOl ITU . 
Because of the E sym (p)5 2 term in the EOS of asymmetric nuclear matter, it is energeti- 
cally more favorable to have a higher isospin asymmetry 5 in the high density region with a 
softer symmetry energy functional E syrn (p). In the supra-normal density region, as shown 
in Fig.[TJ the symmetry energy changes from being soft to stiff when the parameter x varies 
from 1 to -2. Thus the value of (n/p) p > po becomes lower as the parameter x changes from 
1 to -2. It is worth mentioning that the initial value of the quantity {n/p) p > Po is about 1.4 
which is less than the average n/p ratio of 1.56 of the reaction system. This is because 
of the neutron-skins of the colliding nuclei, especially that of the projectile 132 Sn. In the 
neutron-rich nuclei, the n/p ratio on the low-density surface is much higher than that in their 
interior. Also because of the E sym (p)5 2 term in the EOS, the isospin-asymmetry in the low 
density region is much lower than the supra-normal density region as long as the symmetry 
increases with density. In fact, as shown in Fig. 2 of Ref. [54], the isospin-asymmetry of 
the low density region can become much higher than the isospin asymmetry of the reaction 
system. 

It is clearly seen that the dense region can become either neutron-richer or neutron- 
poorer with respect to the initial state depending on the symmetry energy functional E sym {p) 
used. As long as the symmetry energy increases with the increasing density, the isospin 
asymmetry of the supra-normal density region is always lower than the isospin asymmetry 
of the reaction system. Thus, even with radioactive beams, the supra-normal density re- 
gion can not be both dense and neutron-rich simultaneously, unlike the situation in the core 
of neutron stars, unless the symmetry energy starts decreasing at high densities. The high 
density behavior of the symmetry energy is probably among the most uncertain properties 
of dense matter as stressed by ll73l 1741 . Indeed, some predictions show that the symmetry 
energy can decrease with increasing density above certain density and may even finally be- 
comes negative. This extreme behavior was first predicted by some microscopic many-body 
theories, see e.g., Refs. ifTTl 1751 1761 . It has also been shown that the symmetry energy can 
become negative at various high densities within the Hartree-Fock approach using the orig- 
inal Gogny force 11771 . the density-dependent M3Y interaction (7811791 and about 2/3 of the 
87 Skyrme interactions that have been widely used in the literature |[80l . The mechanism 
leading to and physical meaning of a negative symmetry energy are the topics of some very 
recent studies [81 ]. It was found that the spin-isopsin dependence of the three-body force 
and the in-medium properties of the short-range tensor force due to the p-meson exchange 
are the main sources of the very uncertain high-density behavior of the nuclear symmetry 
energy. 

The isospin effects in heavy-ion reactions are determined mainly by the E sym {p)5 2 
term in the EOS. One expects a larger effect if the isospin-asymmetry is higher. Thus, 
ideally, one would like to have situations where both the density and isospin asymmetry 
are sufficiently high simultaneously as in the cores of neutron stars in order to observe the 
strongest effects due to the symmetry energy at supra-normal densities. However, since it 
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Figure 3: (Color online) Pressure as a function of density for symmetric (upper panel) and 
pure neutron (lower panel) matter. The green area in the upper panel is the experimental 
constraint on symmetric matter extracted by Danielewicz, Lacey and Lynch [ 8 ] from ana- 
lyzing the collective flow in relativistic heavy-ion collisions. The corresponding constraint 
on the pressure of pure neutron matter, obtained by combining the flow data and an extrap- 
olation of the symmetry energy functionals constrained below 1.2po (Po = 0.16/m -3 ) by 
the isospin diffusion data, is the shaded black area in the lower panel. Results are taken 
partially from Ref. (SJ. 



is the product of the symmetry energy and the isospin-asymmetry that matters, one can still 
probe the symmetry energy at high densities where the isospin asymmetry is generally low 
with symmetry energy functionals that increase with density. Therefore, even if the high 
density region may not be as neutron-rich as in neutron stars, heavy-ion collisions can still 
be used to probe the symmetry energy at high densities useful for studying properties of 
neutron stars. 

For many astrophysical studies (as those in this chapter), it is more convenient to ex- 
press the EOS in terms of the pressure as a function of density and isospin asymmetry. 
In Fig. [3] we show pressure as a function of density for two extreme cases: symmetric 
(upper panel) and pure neutron matter (lower panel). The green area in the density range 
of po = [2.0 — 4.6] is the experimental constraint on the pressure Pq of symmetric nu- 
clear matter extracted by Danielewicz, Lacey and Lynch from analyzing the collective 
flow data from relativistic heavy-ion collisions (8). The pressure of pure neutron matter 
Ppnm = -Po + p 2 dE sym {p) / dp depends on the density behavior of the nuclear symmetry 
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Table 2: Saturation properties of the nuclear EOSs (for symmetric nuclear matter) shown 
in Fig. [3] 



EOS 


Po 


E. 


K 


m*(p ) 


E S ym(Po) 




(fm~ 3 ) 


(MeV) 


{MeV) 


(MeV/c 2 ) 


(MeV) 


MDI 


0.160 


-16.08 


211.00 


629.08 


31.62 


APR 


0.160 


-16.00 


266.00 


657.25 


32.60 


DBHF+Bonn B 


0.185 


-16.14 


259.04 


610.30 


33.71 



The first column identifies the equation of state. The remaining columns exhibit the following quan- 
tities at the nuclear saturation density: saturation (baryon) density; energy-per-particle; compression 
modulus; nucleon effective mass; symmetry energy. 

energy E sym (p). Since the experimental constraints on the symmetry energy from terres- 
trial laboratory experiments are available mostly for densities less than about 1.2po as indi- 
cated by the green and red squares in the lower panel, which is in contrast to the constraint 
on the symmetric EOS that is only available at much higher densities, the most reliable 
estimate of the EOS of neutron-rich matter can thus be obtained by extrapolating the under- 
lying model EOS for symmetric matter and the symmetry energy in their respective density 
ranges to all densities. Shown by the shaded black area in the lower panel is the resulting 
best estimate of the pressure of high density pure neutron matter based on the predictions 
from the MDI interaction with x = and x = — 1 as the lower and upper bounds on the 
symmetry energy and the flow-constrained symmetric EOS. As one expects and consistent 
with the estimate in Ref . [ 8 ] , the estimated error bars of the high density pure neutron matter 
EOS are much wider than the uncertainty range of the symmetric EOS. For the four interac- 
tions indicated in the figure, their predicted EOSs cannot be distinguished by the estimated 
constraint on the high density pure neutron matter. In addition to the MDI EOS, in Fig. [3] 
we show results by Akmal et al. (821 with the A18 + 5v + UIX* interaction (APR) and 
recent Dirac-Brueckner-Hartree-Fock (DBHF) calculations ll83l with Bonn B One-Boson- 
Exchange (OBE) potential (DBHF+Bonn B) [85 ]. (Older calculations of the DBHF+Bonn 
B EOS can be found in Refs. ll84llTT1l .) The saturation properties of the nuclear equations 
of state applied in this work are summarized in Table 12 



3 Equations determining the structure of neutron stars 

For completeness, in the following we briefly recall the equations determining the structure 
of both static and (rapidly) rotating neutron stars based on the well-established formalism 
in the astrophysical literature. As already mentioned in the introduction, neutron stars are 
objects of extremely compressed matter and therefore proper understanding of their proper- 
ties requires application of both General Relativity and the theories of dense matter. In this 
respect, neutron stars provide a direct link between two of the frontiers of modern physics 
- General Relativity and strong interactions in dense matter |4l. The connection between 
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both branches of physics is provided by Einstein's field equations 




(7) 



(/i, v = 0, 1, 2, 3) which couple the Einstein curvature tensor, G^ u , to the energy-momentum 
tensor, 



of stellar matter. In the above equations P and e denote pressure and mass energy den- 
sity, while R^ h ', g^ v ', and R denote the Ricci tensor, the metric tensor, and the Ricci scalar 
curvature respectively(see e.g., |T]). In Eq. ([U) u M is the unit time-like four-velocity satis- 
fying u^u^ = —1. The tensor T^ u contains the EOS of stellar matter in the form -P(e). 
In general, Einstein's field equations and those of the nuclear many-body problem were to 
be solved simultaneously since the baryons and quarks follow the geodesies of the curved 
space-time whose geometry, determined by the Einstein's field equations, is coupled to 
the total mass energy density of matter [4]. In the case of neutron stars, as for all astro- 
physical situations for which the long-range gravitational force can be separated from the 
short-range strong force, the deviation from flat space-time at the length-scale of the strong 
interactions (~ lfm) is practically zero up to the highest densities achieved in the neu- 
tron star interiors. (This is not to be confused with the global length-scale of neutron stars 
(~ 10km) for which M/R ~ 0.3 depending on the star's mass (in units c = G = 1 so 
that Mq rj 1.475fcm).) In other words, gravity curves space-time only on a macroscopic 
scale but to a very good approximation leaves it flat on a microscopic scale. To achieve an 
appreciable curvature on a microscopic level at which the strong interactions dominate the 
particle dynamics mass densities greater than ~ 10 40 g cm~ 3 would be necessary ||4l [86l. 
Under this circumstances the problem of constructing models of neutron stars separates 
into two distinct tasks. First, the short-range effects of the nuclear forces are described by 
the principles of many-body nuclear physics in a local inertial frame (co-moving proper 
reference frame) in which space-time is flat. Second, the coupling between the long-range 
gravitational force and matter is accounted for by solving the general relativistic equations 
for the gravitational field described by the curvature of space-time, leading to the global 
structure of stellar configurations. 

3.1 Static stars 

In the case of spherically symmetric static (non-rotating) stars the metric has the famous 
Schwarzschild form: 



= (e + P)u»u v + Pg^ 



(8) 



ds 2 



dt 2 + e 2K{r) dr 2 + r 2 {d0 2 + sin 2 0# 2 ) 



(9) 



(c = G = 1) where the metric functions 4>{r) and A(r) are given by: 



e 2A(r) = (1 _ 7(r)) 



(10) 



e 2*(r) 



= e 



2A(r) 



(1 - 7 (r)) r > R s tar 



(11) 
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with 



7(r) 



2m(r) 



2M 



if r < R star 



(12) 



if r > R 



star 



For a static star Einstein's field equations (Eq. (O) reduce then to the familiar Tolman- 
Oppenheimer-Volkoff equation (TOV) ll87l [ 



dP(r) e(r)m(r) 



dr 



1 + 



P(r) 
e(r) 



1 + 



4-7rr 3 p(r) 
m(r) 



1 



2m(r) 



(13) 



where the gravitational mass within a sphere of radius r is determined by 

dm(r) 



dr 



47re(r)r 



The metric function 4>{r) is determined through the following differential equation: 

#(r) __ 1 dP(r) 

dr e(r) + P(r) dr ' 

with the boundary condition at r = R 

1 / 2M\ 
0(r = il) = ~ln(l -— J 



(14) 



(15) 



(16) 



To proceed to the solution of these equations, it is necessary to provide the EOS of 
stellar matter in the form P(e). Starting from some central energy density e c = e(0) at the 
center of the star (r = 0), and with the initial condition m(0) = 0, the above equations 
can be integrated outward until the pressure vanishes, signifying that the stellar edge is 
reached. Some care should be taken at r = since, as seen above, the TOV equation is 
singular there. The point r = R where the pressure vanishes defines the radius of the star 
and M = m(R) = 47r e{r')r' 2 dr' its gravitational mass. 

For a given EOS, there is a unique relationship between the stellar mass and the central 
density e c . Thus, for a particular EOS, there is a unique sequence of stars parameterized by 
the central density (or equivalently the central pressure -P(O)). 



3.2 Rotating stars 

Equations of stellar structure of (rapidly) rotating neutron stars are considerably more com- 
plex than those of spherically symmetric stars (H. These complications arise due to the 
rotational deformations in rotating stars (i.e., flattening at the poles and bulging at the equa- 
tor), which lead to a dependence of the star's metric on the polar coordinate 9. In addition, 
rotation stabilizes the star against gravitational collapse and therefore rotating neuron stars 
are more massive than static ones. A larger mass, however, causes greater curvature of 
space-time. This renders the metric functions frequency-dependent. Finally, the general rel- 
ativistic effect of dragging the local inertial frames implies the occurrence of an additional 
non-diagonal term, g^, in the metric tensor g^ w . This term imposes a self-consistency con- 
dition on the stellar structure equations, since the degree at which the local inertial frames 
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are dragged along by the star, is determined by the initially unknown stellar properties like 
mass and rotational frequency J4J. 

The structure equations of rapidly rotating neutron stars have been computed with the 
-RiVsQcode developed and made available to the public by Nikolaos Stergioulas [46 ]. Here 
we outline briefly the equations solved by the RNS code. The coordinates of the station- 
ary, axial symmetric space-time used to model a (rapidly) rotating neutron star are defined 
through a generalization of Bardeen's metric 11511 : 

ds 2 = -e^ +p dt 2 + e 2a {dr 2 + r 2 d6 2 ) 

+ e 7_p r 2 sin 2 £?(# - ujdt 2 ), (17) 



where the metric potentials 7, p, a, and the angular velocity of the stellar fluid relative to 
the local inertial frame, uj, are functions of the quasi-isotropic radial coordinates, r, and the 
polar angle 9 only. The matter inside a rigidly rotating star is approximated as a perfect 
fluid [51], whose energy momentum tensor is given by Eq. ((8]). The proper velocity of 
matter, v, relative to the local Zero Angular Momentum Observer (ZAMO) [89 ] is defined 
as 

v = rsin(#)(ft-w) e - p(r) (18) 
with £1 = u 3 /u° the angular velocity of a fluid element. The four- velocity is given by 

( 7+ p)/2 

= (1,0,0,0) (19) 



V 



2\ 



In the above equation the function (7 + p) /2 represents the relativistic generalization of the 
Newtonian gravitational potential, while exp[(7 + p) /2] is a time dilation factor between an 
observer moving with angular velocity uj and one at infinity. Substitution of Eq. ( fT9l into 
Einstein's fields equations projected onto the ZAMO reference frame gives three elliptic 
partial differential equations for the metric potentials 7, p, and uj, and two linear ordinary 
differential equations for the metric potential a. Technically, the elliptic differential equa- 
tions for the metric functions are converted into integral equations which are then solved 
iteratively applying Green's function approach lEflllBTl . 

From the relativistic equations of motion, the equations of hydrostatic equilibrium for a 
barotropic fluid may be obtained as [|5T] [ 



h{P) -h p = ) i [uj p + p p - 1 - p- ln(l - v 2 ) + A 2 (uj - n c ) 2 ], (20) 

with h(P) the specific enthalpy, P p the re-scaled pressure, h p the specific enthalpy at the 
pole, 7p and p p the values of the metric potentials at the pole, Q, c = r e f2, and A a rotational 
constant I5T1 . The subscripts p, e, and c label the corresponding quantities at the pole, 
equator and center respectively. The RNS code solves iteratively the integral equations 
for p, 7 and uj, and the ordinary differential equation for the metric function a coupled 
with Eq. (1201 and the equations for hydrostatic equilibrium at the stellar center and equator 
(given h(P c ) and h(P e ) = 0) to obtain p, 7, a, uj, the equatorial coordinate radius r e , 
angular velocity U, energy density e, and pressure P throughout the star. 



'Thanks to Nikolaos Stergioulas the RNS code is available as a public domain program at 
http://www.gravity.phys.uwm.edu/rns/ 
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Figure 4: (Color online) Neutron star masses and radii. Neutron star equatorial radii (upper 
panel) and total gravitational mass (lower panel) versus central energy density e c . Both 
static (solid lines) and Keplerian (broken lines) models are shown. Taken from Ref. [59]. 

4 Constraining global properties of pulsars 

We construct one-parameter 2-D stationary configurations of (rapidly) rotating neutron stars 
employing several nucleonic EOSs and the RNS code. We assume a simple model of stel- 
lar matter of nucleons and light leptons (electrons and muons) in beta-equilibrium. Below 
the baryon density of approximately 0.07/m -3 the equations of state applied here are sup- 
plemented by a crustal EOS, which is more suitable for the low density regime. Namely, 
we apply the EOS by Pethick et al. 1901 for the inner crust and the one by Haensel and 
Pichon ||9~Tll for the outer crust. In the core we assume a continuous functional for the 
EOSs employed in this work. (See Ref. ifTTTl for a detailed description of the extrapolation 
procedure for the DBHF+Bonn B EOS.) 

4.1 Keplerian (and static) sequences 

In what follows we examine the effect of ultra-fast rotation at the Kepler frequency on the 
neutron star gravitational masses and radii. The equilibrium configurations for both static 
and (rapidly) rotating neutron stars are parameterized in terms of the central mass energy 
density, e c = e(0) (or equivalently central pressure, P c = -P(O)). This functional depen- 
dence is shown in Fig. HI where we display the stellar equatorial radius (upper frame) and 
total gravitational mass (lower frame) versus central energy density for the EOSs applied in 
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Figure 5 : (Color online) Mass-radius relation. Both static (solid lines) and Keplerian (bro- 
ken lines) sequences are shown. The 1 — a error bar corresponds to the measurement of the 
mass and radius of EXO 0748-676 gSJ. Taken from Ref. l59l . 




Figure 6: (Color online) Neutron-star polar radius versus central energy density. Both static 
(solid lines) and Keplerian (broken lines) sequences are shown. The labeling of the curves 
is the same as in Fig.|U Taken from Ref. ||59ll . 



this work. Predictions for both static and maximally rotating models are shown. The main 
feature of the mass-density plot is that there exists a maximum value of the gravitational 
mass of a neutron star that a given EOS can support |@). This holds for both static and 
(rapidly) rotating stars. The sequences shown in Fig. [4] terminate at the "maximum mass" 
point. Comparing the results for static and rotating stars, it is seen clearly that the rapid rota- 
tion increases noticeably the mass that can be supported against collapse while lowering the 
central density of the maximum-mass configuration. This is what one should expect, since, 
as already mentioned, rotation stabilizes the star against the gravitational pull providing an 
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Table 3: Maximum-mass static (non-rotating) models. 



EOS 


M max {M Q ) 


R(km) 


e c (xl0 15 <? cm 3 ) 


MDI(x=0) 


1.91 


9.89 


3.02 


APR 


2.19 


9.98 


2.73 


MDI(x=-l) 


1.97 


10.85 


2.57 


DBHF+Bonn B 


2.26 


10.91 


2.37 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities for the static models with maximum gravitational mass: gravitational mass; radius; central 
mass energy density. 

Table 4: Maximum-mass rapidly rotating models at the Kepler frequency v = 



EOS 


M max (M Q ) 


Increase (%) 


e c (xl0 15 g cm 3 ) 






MDI(x=0) 


2.25 


15 


2.59 


1742 


1610 


APR 


2.61 


17 


2.53 


1963 


1699 


MDI(x=-l) 


2.30 


14 


2.21 


1512 


1423 


DBHF+Bonn B 


2.71 


17 


1.94 


1644 


1510 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities for the maximally rotating models with maximum gravitational mass: gravitational mass; 
its percentage increase over the maximum gravitational mass of static models; central mass energy 
density; maximum rotational frequency; Kepler frequency as computed via the empirical relation 
given by Eq. (EB |93l 



extra (centrifugal) repulsion. The rotational effect on the mass-radius relation is illustrated 
in Fig. [5] where the gravitational mass is given as a function of the circumferential radius. 
For rapid rotation at the Kepler frequency, a mass increase up to ~ 17% (Table [4]) is ob- 
tained, depending on the EOS. The equatorial radius increases by several kilometers, while 
the polar radius decreases by several kilometers (see Fig. [6]) leading to an overall oblate 
shape of the rotating star. Table |3] summarizes the properties (masses, radii and central en- 
ergy densities) of the maximum-mass non-rotating neutron star configurations. Our studies 
on the effect of rapid rotation on the upper mass limits for the four EOSs considered in this 
chapter are presented in Table HI In each case the upper mass limit is obtained for a model 
at the mass-shedding limit where v = with central density ~ 15% below that of the 
static model with the largest mass. These findings are consistent with those in Refs. 11921 
and B6l . Table @] also provides an estimate of the upper limiting rotation rate of a neutron 
star. In general, softer EOSs permit larger rotational frequencies since the resulting stellar 
models are more centrally condensed, see e.g., Ref. ll92l . In the last column of Table [4] we 
show the Kepler frequencies computed via the empirical relation 
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Figure 7: (Color online) Mass versus Keplerian (mass-shedding) frequency v\.. Adapted 
fromRef. Wi. 

proposed in Ref. ll93l . The uncertainty of Eq. (|2TT) is ~ 10% (see Ref. ll94l for an improved 
version of the empirical formula). At the time of constructing the above relation Friedman 
et al. [93] did not consider a then unknown class of minimum period EOSs ll96l which 
explains why the numbers in the last column of Table |4] exhibit larger deviation from the 
exact numerical solutions (in column five). This is particularly pronounced for the APR 
EOS for which the approximated Kepler frequency deviates ~ 14% from the exact solution. 
(Note that the APR EOS has the lowest period among the EOSs considered here.) 

Fig.|7]displays the neutron star gravitational mass as a function of the Kepler frequency, 
Vk- The vertical doted line corresponds to the frequency of the fastest neutron star presently 
known (J1748-2446ad 0). The figure suggests that the higher the rotational frequency, the 
larger the mass for the possible stable neutron star configurations. We discuss this further 
in the next subsection. 

4.2 Rotation at various frequencies 

In this subsection we study neutron stars rotating at various (fixed) frequencies. Stability 
with respect to the mass-shedding from equator implies that at a given gravitational mass the 
equatorial radius R eq should be smaller than R™ ax corresponding to the Keplerian limit Q. 
The value of R™ q ax results from the condition that the frequency of a test particle at circular 
equatorial orbit of radius i?™ ax ' just above the equator of the actual rotating star is equal to 
the rotational frequency of the star. As reported by Bejger et al. [2] the relation between 
M and R eq at the "mass-shedding point" is very well approximated by the expression for 
the orbital frequency for a test particle orbiting at r = R eq in the Schwarzschild space-time 
created by a spherical mass. The formula satisfying v^t™' = v, represented by the dotted 
line in Figs. [8l l9land[T0lis given by 




(22) 
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Figure 8: (Color online) Gravitational mass versus circumferential radius for neutron stars 
rotating a.tu = 6A2Hz. See text for details. 
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Figure 9: (Color online) Gravitational mass versus circumferential radius for neutron stars 
rotating at v = 716Hz. 



where v = 6A2Hz in Fig. [8j v = 7WHz in Fig. |9j and v = 1000Hz in Fig. \T0\ respec- 
tively. (Rotational frequencies 6A2Hz and 716Hz are the spinning rates of the two fastest 
pulsars PSR1937+214 and J1748-2446ad 0.) This formula for the Schwarzschild met- 
ric coincides with the one obtained in Newtonian gravity for a point mass M 0. Eq. (l22l 
implies 

/ M \ 1/3 

Rmax = X — nnr- km, (23) 
A \1AM G J ' V ' 

with x = 22.52 for rotational frequency v = 6A2Hz (Fig. U), x = 20.94 for v = 716Hz 
(Fig.©, and X = 16.76 for v = WOOHz (Fig. [TO]). 

In Figs.[8lfT0lwe observe that the range of the allowed masses supported by a given EOS 
for rapidly rotating neutron stars becomes narrower than the one of static configurations. 
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Figure 10: (Color online) Gravitational mass versus circumferential radius for neutron star 
models rotating at v = 1000-ff z. 



This effect becomes stronger with increasing frequency and depends upon the EOS. For 
instance, for models rotating at 1000-ffz (Fig. ITOb for the x = —I EOS the allowed mass 
range is ~ O.2M . As already discussed in the literature [2], this observation could explain 
why such rapidly rotating neutron stars are so rare - their allowed masses fall within a very 
narrow range. 



5 Neutron star moment of inertia 



In this section we turn our attention to the moment of inertia of neutron stars. Such stud- 
ies are important and timely as they are related to astrophysical observations in the near 
future. In particular, the moment of inertia of pulsar A in the extremely relativistic neu- 
tron star binary PSR J0737-3039 11971 may be determined in a few years through detailed 
measurements of the periastron advance [98]. 

Employing the EOSs described briefly in section 2, we compute the neutron star mo- 
ment of inertia with the RNS code. It solves the hydrostatic and Einstein's field equations 
for mass distributions rotating rigidly under the assumption of stationary and axial symme- 
try about the rotational axis, and reflectional symmetry about the equatorial plane. RNS 
calculates the angular momentum J as 11511 

J = J T^tfadV, (24) 
where T^ v is the energy-momentum tensor of stellar matter 

T» v = (e + P)u"u v + iV, (25) 

is the Killing vector in azimuthal direction reflecting axial symmetry, and dV = 
y/—gd 3 x is a proper 3-volume element (g = det(g a p) is the determinant of the 3-metric). 
In Eq. (l25l) P is the pressure, e is the mass-energy density, and n M is the unit time-like 
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four-velocity satisfying u^u^ = —1. For axial-symmetric stars it takes the form = 
ti'(l, 0, 0, f2), where Q is the star's angular velocity. Under this condition Eq. (1241) reduces 
to 

J = J(e + P)u t (g H u' t ' + g^ t u t )^d 3 x (26) 

It should be noted that the moment of inertia cannot be calculated directly as an integral 
quantity over the source of gravitational field 1511 . In addition, there exists no unique 
generalization of the Newtonian definition of the moment of inertia in General Relativity 
and therefore I = J/Q is a natural choice for calculating this important quantity. 

For rotational frequencies much lower than the Kepler frequency (the highest possible 
rotational rate supported by a given EOS), i.e. u/uf, << 1 (y = J7/(27r)), the deviations 
from spherical symmetry are very small, so that the moment of inertia can be approximated 
from spherical stellar models. In what follows we review briefly this slow-rotation approxi- 
mation, see e.g. Ref. Il39l . In the slow-rotational limit the metric can be written in spherical 
coordinates as (in geometrized units G = c = 1) 

ds 2 = -e 2 ^ r) dt 2 + ( 1 - dr 2 -2ior 2 sin 2 9dtd0 + r 2 (d8 2 + sin 2 9d(j) 2 ) (27) 



In the above equation m(r) is the total gravitational mass within radius r satisfying the 
usual equation 

Amir) n 

4vre(r)r 2 (28) 



dr 

and w(r) = (dcj)/ dt) zamo is the Lense-Thirring angular velocity of a zero-angular-momentum 
observer (ZAMO). Up to first order in co all metric functions remain spherically symmetric 
and depend only on r (99). In the stellar interior the Einstein's field equations reduce to 



#(r) 

— - — = m(r) 
dr 



1 + 



47rr 3 P(r) |" 1 2m(r)~ l ' 
m{r) 



1 



(r < R star ) (29) 



and 



with u = Q — lo the dragging angular velocity (the angular velocity of the star relative to a 
local inertial frame rotating at u) and 

AMoy'v^, (31) 



Outside the star the metric functions become 

IM 
r J 



e 2 f=(l-—) (r>R star ) (32) 



and 

2J 

W = (r > R s tar), (33) 



where M = m(r = R) = An e{r')r' 2 dr' is the total gravitational mass and R is the 
stellar radius defined as the radius at which the pressure drops to zero (P(r = R) = 0). At 
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the star's surface the interior and exterior solutions are matched by satisfying the appropriate 
boundary conditions 

(34) 



~ Rfduj\ 

q(r) = n — — 

v ; 3 \dr r =R 



and 



(r) 



1 



In 1 



2M\ 



(35) 



The moment of inertia I = J/U then can be computed from Eq. (|26T ). With SI = u^/u 1 and 
retaining only first order terms in to and Q, the moment of inertia reads 



o Jo 



1 



2m(r) 
r 



to 

n 



r dr 



(36) 



This slow-rotation approximation for the neutron-star moment of inertia neglects devia- 
tions from spherical symmetry and is independent of the angular velocity Q, [99]. For 
neutron stars with masses greater than IM & Lattimer and Schutz HI 001 found that, for 
slow-rotations, the moments of inertia computed through the above formalism (Eq. (l36l » 
can be approximated very well by the following empirical relation: 



(0.237 ±0.008)M,R 2 



l + 4.2 MATO +90f M ^V' 

m q r \m q rJ 



(37) 



The above equation is shown H100I to hold for a wide class of EOSs except for ones with 
appreciable degree of softening, usually indicated by achieving a maximum mass of ~ 
1.6M Q or less. Since none of the EOSs employed in this paper exhibit such pronounced 
softening, Eq. d37l ) is a good approximation for the momenta of inertia of slowly rotating 
stars. 



5.1 Slow rotation 

If the rotational frequency is much smaller than the Kepler frequency, the deviations from 
spherical symmetry are negligible and the moment of inertia can be calculated applying 
the slow-rotation approximation discussed briefly above. For this case Lattimer and Schutz 
[100] showed that the moment of inertia can be very well approximated by Eq. (I37T ). In 
Fig. [EU we display the moment of inertia as a function of stellar mass for slowly rotat- 
ing neutron stars as computed with the empirical relation Eq. (l37l) . As shown in Fig. [5J 
above ~ l.OM the neutron star radius remains approximately constant before reaching 
the maximum mass supported by a given EOS. The moment of inertia (I ~ MR 2 ) thus 
increases almost linearly with stellar mass for all models. Right before the maximum mass 
is achieved, the neutron star radius starts to decrease (Fig. [5]), which causes the sharp drop 
in the moment of inertia observed in Fig. [TTJ Since I is proportional to the mass and the 
square of the radius, it is more sensitive to the density dependence of the nuclear symmetry 
energy, which determines the neutron star radius. Here we recall that the x = — 1 EOS has 
much stiffer symmetry energy (with respect to the one of the x = EOS), which results 
in neutron star models with larger radii and, in turn, momenta of inertia. For instance, for 
a "canonical" neutron star (M = 1.4M©), the difference in the moment of inertia is more 
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Figure 11: (Color online) Total moment of inertia of neutron stars estimated with Eq. (I37T ). 
Taken from Ref. (601. 
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Figure 12: (Color online) The moment of inertia scaled by M 3 / 2 as a function of the stellar 
mass M. The shaded band illustrates a 10% error of hypothetical I/M 3 / 2 measurement of 
50 km 2 M~ 1/2 . The error bar shows the specific case in which the mass is 1.34M (after 
Lattimer and Schutz ifTOOln . Taken from Ref. Il60l . 

than 30% with the x = and the x = — 1 EOSs. In Fig. [12] we take another view of the 
moment of inertia where / is scaled by M 3 / 2 as a function of the stellar mass (after [ 100]). 

The discovery of the extremely relativistic binary pulsar PSR J0737-3039A,B provides 
an unprecedented opportunity to test General Relativity and physics of pulsars [97]. Lat- 
timer and Schutz II 1001 estimated that the moment of inertia of the A component of the 
system should be measurable with an accuracy of about 10%. Given that the masses of 
both stars are already accurately determined by observations, a measurement of the mo- 




DBHF+Bonn B 




24 



Plamen G. Krastev and Bao-An Li 



Table 5: Numerical results for PSR J0737-3039A (M A = 1.338M©, v A = 44.05.Hz). 



EOS 


e c (xl0 14 g cm 3 ) 


R eq (km) 


I(xW A5 gcm 2 ) 


I LS (xl0 45 g cm 2 ) 


MDI(x=-l) 


7.04 


13.75 (13.64) 


1.63 


1.67 


DBHF+Bonn B 


7.34 


12.56(12.47) 


1.57 


1.43 


MDI(x=0) 


9.85 


12.00(11.90) 


1.30 


1.34 


APR 


9.58 


11.60(11.52) 


1.25 


1.26 



The first column identifies the equation of state. The remaining columns exhibit the following 
quantities: central mass-energy density, equatorial radius (the numbers in the parenthesis are the 
radii of the spherical models; the deviations from sphericity due to rotation are ~ 1%), total moment 
of inertia, total moment of inertia I LS as computed with Eq. ( TTTt . 



ment of inertia of even one neutron star could have enormous importance for the neutron 
star physics [100]. (The significance of such a measurement is illustrated in Fig. [T2] As 
pointed by H1001 . it is clear that very few EOSs would survive these constraints.) Thus, the- 
oretical predictions of the moment of inertia are very timely. Calculations of the moment of 
inertia of pulsar A (M A = 1.338M Q , v A = 44.05Hz) have been reported by Morrison et 
al. ll99l and Bejger et al. ll98l . In Table[5]we show the moment of inertia (and other selected 
quantities) of PSR J0737-3039A computed with the RNS code using the EOSs employed 
in this study. Our results with the APR EOS are in very good agreement with those by ll99l 
(jAPR = 1 24 x io 4 5 5 cm 2 ) and EH (I APR = 1.23 x 10 45 g cm 2 ). In the last column 
of Table [5] we also include results computed with the empirical relation (Eq. (l37l)). From a 
comparison with the results from the exact numerical calculation we conclude that Eq. (l37l ) 
is an excellent approximation for the moment of inertia of slowly-rotating neutron stars. 
(The average uncertainty of Eq. (l37l) is ~ 2%, except for the DBHF+BonnB EOS for which 
it is ~ 8%.) Our results (with the MDI EOS) allowed us to constrain the moment of inertia 
of pulsar A to be in the range / = (1.30 - 1.63) x 10 45 (g cm 2 ). 

5.2 Rapid rotation 

In this subsection we turn our attention to the moment of inertia of rapidly rotating neutron 
stars. In Fig. [13] we show the moment of inertia as a function of stellar mass for neutron 
star models spinning at the mass-shedding (Kepler) frequency. The numerical calculation is 
performed with the RNS code. We observe that the moments of inertia of rapidly rotating 
neutron stars are significantly larger than those of slowly rotating models (for a fixed mass). 
This is easily understood in terms of the increased (equatorial) radius (Fig. [5]). 

We also compute the moments of inertia of pulsars spinning at 642Hz, 715Hz, and 
1000Hz (642Hz [3 ] and 716Hz [6] are the rotational frequencies of the fastest pulsars as 
of today). The numerical results are presented in Fig. |T4j As demonstrated by Bejger et 
al. O and most recently by Krastev et al. ||59l , the range of the allowed masses supported 
by a given EOS for rapidly rotating neutron stars becomes narrower than the one for static 
configurations. The effect becomes stronger with increasing frequency and depends upon 
the EOS. This is also illustrated in Fig. [Q] particularly in panel (c). Additionally, the 
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Figure 13: (Color online) Total moment of inertia for Keplerian models. The neutron star 
sequences are computed with the RNS code. Taken from Ref . ll60l . 



4-10 



45 



45 



3-10 



- 2-1 45 



1-10 
4- 10 



: v = 642Hz 


— APR 


— x = 




x = -1 




- DBHF+Bonn B - 











v = 716Hz 



3-10 



,45 ; 



- 2-1 45 



1-10 



4-10 




v = 1000Hz 



3-10 



,45- 



- 2-10 



,45 



1-10 



45 [ v ' - 

1.2 1.4 1.6 1.8 2.0 2.2 2.4 

M/M e 



Figure 14: (Color online) Total moment of inertia as a function of stellar mass for models 
rotating at 642Hz (a), 716Hz (b), and 1000Hz (c). Data is partially taken from Ref. [60]. 
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Figure 15: (Color online) Total moment of inertia as a function of rotational frequency for 
stellar models with mass M = 1.4M Q . Taken from Ref. [60]. 



moment of inertia shows increase with rotational frequency at a rate dependent upon the 
details of the EOS. This is best seen in Fig. Q3] where we display the moment of inertia as a 
function of the rotational frequency for stellar models with a fixed mass (M = 1.4M ). 

The neutron star sequences shown in Fig. [14] are terminated at the mass-shedding fre- 
quency. At the lowest frequencies the moment of inertia remains roughly constant for all 
EOSs (which justifies the application of the slow-rotation approximation and Eq. (l37l)). As 
the stellar models approach the Kepler frequency, the moment of inertia exhibits a sharp 
rise. This is attributed to the large increase of the circumferential radius as the star ap- 
proaches the "mass-shedding point". As pointed by Friedman et al. [92], properties of 
rapidly rotating neutron stars display greater deviations from those of spherically symmet- 
ric (static) stars for models computed with stiffer EOSs. This is because such models are 
less centrally condensed and gravitationally bound. This also explains why the momenta 
of inertia of rapidly rotating neuron star configurations from the x = — 1 EOS show the 
greatest deviation from those of static models. 

5.3 Fractional moment of inertia of the neutron star crust 

As it was discussed extensively by Lattimer and Prakash lf35l (and others), the neutron star 
crust thickness might be measurable from observations of pulsar glitches, the occasional 
disrupts of the otherwise extremely regular pulsation from magnetized, rotating neutron 
stars. The canonical model of Link et al. H 1 1 II suggests that glitches are due to the angular 
momentum transfer from superfluid neutrons to normal matter in the neutron star crust, the 
region of the star containing nuclei and nucleons that have dripped out of nuclei. This region 
is bound by the neutron drip density at which nuclei merge into uniform nucleonic matter. 
Link et al. H 1 1 II concluded from the observations of the Vela pulsar that at least 1 .4% of the 
total moment of inertia resides in the crust of the Vela pulsar. For slowly rotating neutron 
stars, applying several realistic hadronic EOSs that permit maximum masses of at least 
~ 1.6M© Lattimer and Prakash OBI found that the fractional moment of inertia, AI/I, 
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Table 6: Transition densities and pressures for the EOSs used in this paper. 



EOS 


MDI(x=0) 


MDI(x=-l) 


APR 


DBHF+Bonn B 


Pt(fm~ 3 ) 


0.091 (0.095) 


0.093 (0.092) 


0.087 


0.100 


P t (MeV fm- 3 ) 


0.645 


0.982 


0.513 


0.393 



The first row identifies the equation of state. The remaining rows exhibit the following quantities: 
transition density, transition pressure. The numbers in the parenthesis are the transition densities 
calculated by Kubis lfT02l . 



can be expressed approximately as 
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In the above equation AI is the moment of inertia of the neutron star crust, I is the total 
moment of inertia, /3 = GM/Rc 2 is the compactness parameter, is the average nucleon 
mass, p t is the transition density at the crust-core boundary, and P t is the transition pressure. 
The determination of the transition density itself is a very complicated problem. Different 
approaches often give quite different results. Similar to determining the critical density for 
the spinodal decomposition for the liquid-gas phase transition in nuclear matter, for uniform 
rape-matter, Lattimer and Prakash ll35l and more recently Kubis 1110211 have evaluated the 
crust transition density by investigating when the incompressibility of rape-matter becomes 
negative, i.e 
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Figure 16: (Color online) The incompressibility, K^, as a function of baryon density p. 
Taken from Ref. flU. 
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Figure 17: (Color online) The fractional moment of inertia of the neutron star crust as 
a function of the neutron star mass (left panel) and radius (right panel) estimated with 
Eq. (l39l >. The constraint from the glitches of the Vela pulsar is also shown. Taken from 
Ref. Il60l . 



(see Fig.[T6l) where Eq(p) is the EOS of symmetric nuclear matter, E sym is the nuclear sym- 
metry energy, and 5 = (p n —Pp) / (Pn+Pp) is the asymmetry parameter. Using this approach 
and the MDI interaction, Kubis H102II found the transition density of 0.119, 0.092, 0.095 and 
0.160/m -3 for the x parameter of 1, 0, —1 and —2, respectively. Similarly, we have cal- 
culated the transition densities and pressures for the EOSs employed in this work. Our 
results are summarized in Table [6] We find good agreement between our results and those 
by Kubis [102] with the MDI interaction. It is interesting to notice that the transition den- 
sities predicted by all EOSs are in the same density range explored by heavy-ion reactions 
at intermediate energies. The MDI interaction with x = and x = — 1 constrained by 
the available data on isospin diffusion in heavy-ion reaction at intermediate energies thus 
limits the transition density rather tightly in the range of Pt = [0.091 - 0.093] (fm- 3 ). 
It is worth noticing, however, that the transition density is estimated here by using the 
parabolic approximation of the EOS for isospin asymmetric nuclear matter. Relaxing this 
approximation, using both dynamical and thermodynamic al approaches, the transition den- 
sity and pressure have been studied in great detail in Ref. [56]. A somewhat lower values 
of p t = [0.040 — 0.065] (fm~ 3 ) were obtained. Nevertheless, to be consistent with the rest 
of the study in this review, we stick to the parabolic approximation of the EOS in studying 
the AI/I. 

The fractional momenta of inertia AI/I of the neutron star crusts are shown in Fig.fTTI 
as computed through Eq. (l38l > with the parameters listed in Table |6] It is seen that the 
condition AI/I > 0.014 extracted from studying the glitches of the Vela pulsar does put a 
strict lower limit on the radius for a given EOS. It also limits the maximum mass to be less 
than about 2M for all of the EOSs considered. Similar to the total momenta of inertia the 
ratio AI/I changes more sensitively with the radius as the EOS is varied. 
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Figure 18: (Color online) Density (upper panel) and proton fraction (lower panel) versus 
rotational frequency for fixed neutron star mass. Taken from Ref. [59]. 

6 Rotation and proton fraction 

Finally, we study the effect of (fast) rotation on the proton fraction in the neutron star core. 
In Fig. [18] we show the central baryon density (upper frame) and central proton fraction 
(lower frame) as a function of the rotational frequency for fixed-mass models. Predictions 
from both x = and x = — 1 EOSs are shown. We observe that central density decreases 
with increasing frequency. This reduction is more pronounced in heavier neutron stars. 
Most importantly, we also observe decrease in the proton fraction Y£ in the star's core. 
We recall that large proton fraction (above ~ 0.14 for npe^i-stars) leads to fast cooling 
of neutron stars through direct Urea reactions. Our results demonstrate that depending on 
the stellar mass and rotational frequency, the central proton fraction could, in principle, 
drop below the threshold for the direct nucleonic Urea channel and thus making the fast 
cooling in rotating neutron stars impossible. The masses of the models shown in Fig.[T8lare 
chosen so that the proton fraction in stellar core is just above the direct Urea limit for the 
static configurations, see Fig. [19] upper and lower frames. The stellar sequences in Fig. IT81 
are terminated at the Kepler (mass-shedding) frequency. In both cases the central proton 
fraction drops below the direct Urea limit at frequencies lower than that of PSR J1748- 
244ad [5]. This implies that the fast cooling can be effectively blocked in millisecond 
pulsars depending on the exact mass and spin rate. It might also explain why heavy neutron 
stars (could) exhibit slow instead of fast cooling. For instance, with the x = EOS (with 
softer symmetry energy) for a neutron star of mass approximately 1.9M Q , the Direct Urea 
channel closes at u A70Hz. On the other hand, with the x=-l EOS (with stiffer symmetry 
energy) the direct Urea channel can close only for low mass neutron stars, in fact only for 
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Figure 19: (Color online) Neutron star mass (upper panel), nuclear symmetry energy (mid- 
dle panel) and proton fraction (lower panel) versus baryon density. Predictions from the 
MDI (x=-l,0) are shown. Taken from Ref. ll59l . 



masses well below the canonical mass of 1.4M©. This is due to the much stiffer symmetry 
energy (see Fig.[l9]middle frame) because of which the direct Urea threshold (Fig.[T9llower 
frame) is reached at much lower densities and stellar masses (Fig. [19] upper frame). 

Before closing the discussion in this section a few comments are in order. In the present 
study we do not consider "exotic" states of matter in neutron stars. On the other hand, 
due to the rapid rise of the baryon chemical potentials several other species of particles, 
such as strange hyperons A and E~, are expected to appear once their mass thresholds 
are reached 11031 . The appearance of hyperonic degrees of freedom lowers the energy- 
per-particle in the stellar medium and causes more centrally condensed configurations with 
lower masses and radii. Additionally, hyperons help the condition for nucleonic direct Urea 
process to be satisfied at lower densities due to the increased proton fraction [104], and 
depending on their exact concentrations could potentially contribute to the fast cooling of 
the star through hyperonic direct Urea processes 111051 . These considerations would alter 
the balance between the curves in Fig. [19] and ultimately the results displaced in Fig. [18] 
in favor of the direct Urea process, i.e. smaller masses and higher frequencies would be 
necessary to close the fast cooling channel. This is due to the fact that the overall impact of 
(rapid) rotation on the neutron star structure is smaller for more centrally condensed models 
resulting from "softer" EOSs ||92ll . Therefore, for such models there is smaller deviation 
from properties and structure of static configurations. In addition, at even higher densities 
matter is expected to undergo a transition to quark-gluon plasma Bl 11031 . which favors a 
fast cooling through enhanced nucleonic direct Urea and quark direct Urea processes, see 
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e.g., Ref. una. 

7 Constraining gravitational waves from elliptically deformed 
pulsars 

Gravitational waves are tiny disturbances in space-time and are a fundamental, although 
not yet directly confirmed, prediction of General Relativity. They can be triggered in cat- 
aclysmic events involving (compact) stars and/or black holes. They could even have been 
produced during the very early Universe, well before any stars had been formed, merely 
as a consequence of the dynamics and expansion of the Universe. Because gravity interact 
extremely weakly with matter, gravitational waves would carry a genuine picture of their 
sources and thus provide undisturbed information that no other messenger can deliver 111061 . 
Gravitational wave astrophysics would open an entirely new non-electromagnetic window 
making it possible to probe physics that is hidden or dark to current electromagnetic obser- 
vations ifTfJTI . 

(Rapidly) rotating neutron stars could be one of the major candidates for sources of con- 
tinuous gravitational waves in the frequency bandwidth of the LIGO [ 108 ] and VIRGO (e.g. 
Ref. Cl09]) laser interferometric detectors. It is well known that a rotating object self-bound 
by gravity and which is perfectly symmetric about the axis of rotation does not emit grav- 
itational waves. In order to generate gravitational radiation over extended period of time, 
a rotating neutron star must have some kind of long-living axial asymmetry [ 1 10]. Several 
mechanisms leading to such an asymmetry have been studied in the literature: (1) Since the 
neutron star crust is solid, its shape might not be necessarily symmetric, as it would be for 
a fluid, with asymmetries supported by anisotropic stress built up during the crystallization 
period of the crust Mill . (2) Additionally, due to its violent formation (supernova) or due to 
its environment (accretion disc), the rotational axis may not coincide with a principal axis 
of the moment of inertia of the neutron star which make the star precess H112L Even if the 
star remains perfectly symmetric about the rotational axis, since it precesses, it emits gravi- 
tational waves nil 21 1 1 1 3ll . (3) Also, the extreme magnetic fields presented in a neutron star 
cause magnetic pressure (Lorenz forces exerted on conducting matter) which can distort the 
star if the magnetic axis is not aligned with the axis of rotation 11141 . which is widely sup- 
posed to occur in order to explain the pulsar phenomenon. Several other mechanisms exist 
that can produce gravitational waves from neutron stars. For instance, accretion of matter on 
a neutron star can drive it into a non-axisymmetric configuration and power steady radiation 
with a considerable amplitude 11151 . This mechanism applies to a certain class of neutron 
stars, including accreting stars in binary systems that have been spun up to the first insta- 
bility point of the so-called Chadrasekhar-Friedman-Schutz (CFS) instability 11161 . Also, 
Andersson 111171 suggested a similar instability in r-modes of (rapidly) rotating relativistic 
stars. It has been shown that the effectiveness of these instabilities depends on the viscosity 
of stellar matter which in turn is determined by the star's temperature. 

Gravitational waves are characterized by tiny dimensionless strain amplitude, which 
depends on the degree to which the neutron star is deformed from axial symmetry which, 
in turn, depends upon the equation of state (EOS) of neutron-rich stellar matter. As already 
mentioned, presently the EOS of matter under extreme conditions (densities, pressures and 
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isospin asymmetries) is still rather uncertain, mainly due to the poorly known density de- 
pendence of the nuclear symmetry energy, E sym (p), e.g. |9]. Applying several nucleonic 
EOSs (discussed in section 2), we calculate the gravitational wave strain amplitude for se- 
lected neutron star configurations. Particular attention is paid to predictions with the EOS 
partially constrained by the nuclear laboratory data. These results set an upper limit on the 
strain amplitude of gravitational radiation expected from rotating neutron stars. 

The pulsar population is such that most have spin frequencies that fall below the sensi- 
tivity band of current detectors. In the future, the low-frequency sensitivity of VIRGO [118] 
and Advanced LIGO 111 191 should allow studies of a significantly larger sample of pulsars. 
Moreover, LISA (the Laser Interferometric Space Antenna) is currently being jointly de- 
signed by NASA in the United States and ESA (the European Space Agency), and will be 
launched into orbit in the near future providing an unprecedented instrument for gravita- 
tional waves search and detection 11071 . The discussion in this section is (mainly) based on 
Ref. ED. 



7.1 Formalism relating the EOS of neutron-rich matter to the strength of 
gravitational waves from slowly rotating neutron stars 

In the following we review briefly the formalism used to calculate the gravitational wave 
strain amplitude from slowly rotating neutron stars. A spinning neutron star is expected 
to emit GWs if it is not perfectly symmetric about its rotational axis. As already men- 
tioned, non-axial asymmetries can be achieved through several mechanisms such as elas- 
tic deformations of the solid crust or core or distortion of the whole star by extremely 
strong misaligned magnetic fields. Such processes generally result in a triaxial neutron star 
configuration [108] which, in the quadrupole approximation and with rotation and angu- 
lar momentum axes aligned, would cause gravitational waves at twice the star's rotational 
frequency 1110811 . These waves have characteristic strain amplitude at the Earth's vicinity 
(assuming an optimal orientation of the rotation axis with respect to the observer) of II 1201 

Wir 2 GeI zz v 2 

ho = 4 — , (40) 

c 4 r 

where v is the neutron star rotational frequency, I zz its principal moment of inertia, e = 
(I xx — I yy )/Izz its equatorial ellipticity, and r its distance to Earth. The ellipticity is related 
to the neutron star maximum quadrupole moment (with m = 2) via 11 1 2 1 II 
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where for slowly rotating (and static) neutron stars $22 can be written as B1211 



(41) 



= 2A y 10389 °»' UcrJ Uo^J l"ir ) (42) 

In the above expression a is the breaking strain of the neutron star crust which is rather 
uncertain at present time. Earlier studies have estimated a to be in the range a = [10 -5 — 
10 -2 ] 1112011 . More recently, using molecular dynamics simulations, it was estimated to be 
about . 1 which is considerably larger than the previous findings 111221 . In our work ||6T1 , we 
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have used a = 10 -2 , which, compared to the latest estimate H122H . is rather conservative. 
Nevertheless, our results are simply amplified by a factor of 100 if we apply the estimate by 
Horowitz et al. HI 2211 in our calculations. From Eqs. (1401 and (|4TT ) it is clear that Iiq does not 
depend on the moment of inertia I zz , and that the total dependence upon the EOS is carried 
by the quadrupole moment $22- Thus Eq. (1401 can be rewritten as 

h = X ^^, (43) 
r 

with x = \/20457r 5 /15G'/c 4 . Eq. (l43l establishes the link between the gravitational 
wave strain amplitude and the EOS of neutron-rich matter-the EOS is the main ingredi- 
ent for determining the neutron star properties, including its quadrupole moment. In a 
recent work ll60l we have calculated the neutron star moment of inertia of both static and 
(rapidly) rotating neutron stars. For slowly rotating neutron stars Lattimer and Schutz 111001 
derived an empirical relation (Eq. (1371) ) for I zz which is shown to hold for a wide class of 
equations of state which do not exhibit considerable softening and for neutron star models 
with masses above 1M Q [100]. Using Eq. (137T ) to calculate the neutron star moment of 
inertia and Eq. (142b the corresponding quadrupole moment, the ellipticity e can be readily 
computed (via Eq. (I4TT) ). Since the global properties of spinning neutron stars (in partic- 
ular the moment of inertia) remain approximately constant for rotating configurations at 
frequencies up to ~ 300 Hz [60], the above formalism can be readily employed to estimate 
the gravitational wave strain amplitude, provided one knows the exact rotational frequency 
and distance to Earth, and that the frequency is relatively low (below ~ 300Hz). These 
estimates are then to be compared with the current upper limits for the sensitivity of the 
laser interferometric observatories (e.g. LIGO). 



7.2 Gravitational waves from slowly rotating neutron stars 

We calculate the gravitational wave strain amplitude Iiq for several selected pulsars with 
rotational frequencies ~ 200Hz and relatively close to Earth, employing several nucleonic 
equations of state. We assume a simple model of stellar matter of nucleons and light leptons 
(electrons and muons) in beta-equilibrium. The equations of state that we employ here have 
been discussed in section 2 and shown in Fig. [3] 

Fig. [20] displays the neutron star quadrupole moment (left panel) and ellipticity (right 
panel). The quadrupole moment is calculated through Eq. (|42l and the ellipticity through 
Eq. (|4TT ). Note that Eq. (l42l is valid only for slowly rotating neutron star models. We 
notice that $22 decreases with increasing stellar mass for all EOSs considered in this study. 
The rate of this decrease depends upon the EOS and is largest for the x = — 1 EOS. This 
behavior is easily understood in terms of the increased central density with stellar mass 
- more massive stars are more compact and, since the quadrupole moment is a measure 
of the star's deformation (see Eq. (PHI), they are also less deformed with respect to less 
centrally condensed models. Moreover, it is well known that the mass is mainly determined 
by the symmetric part of the EOS while the radius of a neutron star is strongly affected 
by the density slope of the symmetry energy. More quantitatively, an EOS with a stiffer 
symmetry energy, such as the x = — 1 EOS, results in less compact stellar models, and 
hence more deformed pulsars. Here we recall specifically that the x = — 1 EOS yields 
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Figure 20: (Color online) Neutron star quadrupole moment (left panel) and ellipticity (right 
panel). Taken from Ref. lloTll . 



neutron star configurations with larger radii than those of models form the rest of the EOSs 
considered in this study (e.g. see Fig. [5]). These results are consistent with previous findings 
which suggest that more compact neutron star models are less altered by rotation, e.g. see 
Ref. [92]. Consequently, it is reasonable also to expect such configurations to be more 
"resistant" to any kind of deformation. The neutron star ellipticity, e, is shown as a function 
of the neutron star mass (Fig. [20j right panel). Since e is proportional to the quadrupole 
moment $22 (scaled by the moment of inertia I zz ), it decreases with increasing stellar mass. 
The results shown in the right panel of Fig. [20] are consistent with the maximum ellipticity 
tmax ~ 2.4 x 10~ 6 corresponding to the largest crust "mountain" one could expect on 
a neutron star II 1201 11231 . (The estimate of e ma x in Refs. 111201 [123 1 has been obtained 
assuming breaking strain of the crust a = 10~ 2 , as we have assumed in the present paper 
in calculating the neutron star quadrupole moment.) 

In Fig.[2T]we display the GW strain amplitude, ho, as a function of stellar mass. Predic- 
tions are shown for three selected millisecond pulsars, which are relatively close to Earth 
(r < OAkpc), and have rotational frequencies below 300Hz so that the corresponding mo- 
menta of inertia and quadrupole moments can be computed approximately via Eqs. (I37T ) and 
(l42l respectively. The properties of these pulsars (of interest to this study) are summarized 
in Table [7] The error bars in Fig. [21] between the x = and x = — 1 EOSs provide a 
constraint on the maximal strain amplitude of the gravitational waves emitted by the mil- 
lisecond pulsars considered here. The specific case shown in the figure is for neutron star 
models of 1.4M . Depending on the exact rotational frequency, distance to detector, and 
details of the EOS, the maximal ho is in the range ~ [0.4—1.5] x 10~ 24 . These estimates do 
not take into account the uncertainties in the distance measurements. They also should be 
regarded as upper limits since the quadrupole moment (Eq. (142b ) has been calculated with 
a = 10 -2 (where a can go as low as 10~ 5 ). Here we recall that the mass of PSR J0437- 
4715 (Fig. EDright panel) is 1.3±O.2M HT25\ . (Another mass constraint, 1.58±O.18M , 
was given previously by van Straten et al. 11261 .) The results shown in Fig. [21] suggest that 
the GW strain amplitude depends on the EOS of stellar matter, where this dependence is 
stronger for lighter neutron star models. In addition, it is also greater for stellar configura- 



Nuclear limits on properties of pulsars and gravitational waves 



35 



3.0i 

: v = 202.79Hz, r = 0.25kpc 

2.5 N. PSRJ21 24-3358 




v = 193.72Hz, r = 0.35kpc 


PSRJ1 024-071 9 : 















v = 173.91Hz, r = 0.18kpc 


PSR J0437-4715 















1.0 1.2 1.4 1.6 1.8 2.0 2.2 1.0 1.2 1.4 1.6 1.8 2.0 2.2 1.0 1.2 1.4 1.6 1.8 2.0 2.2 



M/IVL 



M/lvL 



M/M, 



Figure 21: (Color online) Gravitational- wave strain amplitude as a function of the neutron 
star mass. The error bars between the x = and x = — 1 EOSs provide a limit on the strain 
amplitude of the gravitational waves to be expected from these neutron stars, and show a 
specific case for stellar models of IAM & . Taken form Ref. 116X11 . 



Table 7: Properties of the pulsars considered in this study. 



Pulsar 


v(Hz) 


M(M ) 


r(kpc) 


Reference 


PSR J2124-3358 


202.79 




0.25 


11241 


PSRJ1024-0719 


193.72 




0.35 


111241 


PSR J0437-4715 


173.91 


1.3 ±0.2 


0.18 


1125111261 



The first column identifies the pulsar. The remaining columns exhibit the following quantities: 
rotational frequency; mass (if known); distance to Earth; corresponding references. Notice that 
only the mass of PSR J0437-4715 is known from orbital dynamics 1 126, 1251 ( as the pulsar has a 
low-mass white dwarf companion). The masses of PSRs J2124-3358 and J1024-0719 are presently 
unknown as they are both isolated neuron stars 01241 . 



tions computed with stiffer EOS. As explained, such models are less compact and thus less 
gravitationally bound. As a result, they could be more easily deformed by rotation or/and 
other deformation driving mechanisms and phenomena, and therefore are expected to emit 
stronger gravitational radiation (Eq. (l40l). 

In Fig. [22] we take another view of the results shown in Fig. [21] We display the maxi- 
mal GW strain amplitude as a function of the GW frequency and compare our predictions 
with the best current detection limit of LIGO. The specific case shown is for neutron star 
models with mass 1.4M Q computed with the x = and x = — 1 EOSs. Since these EOSs 
are constrained by the available nuclear laboratory data they provide a limit on the possible 
neutron star configurations and thus gravitational emission from them. The results shown 
in Fig. [22] would suggest that presently the gravitational radiation from the three selected 
pulsars should be within the detection capabilities of LIGO. The fact that such a detection 
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Figure 22: (Color online) Gravitational wave strain amplitude as a function of the gravita- 
tional wave frequency. The characters denote the strain amplitude of the GWs expected to 
be emitted from spinning neutron stars (v < 300Hz) with mass 1.4M . Solid line denotes 
the current upper limit of the LIGO sensitivity. Adapted from Ref. H108B . 

has not been made yet deserves a few comments at this point. First, as we mentioned, in 
the present calculation we assume breaking strain of the neutron star crust a = 10~ 2 which 
might be too optimistic. As pointed by Haskell et al. II 1201 . we are still away from testing 
more conservative models and if the true value of a lies in the low end of its range, we 
would be still far away from a direct detection of a gravitational wave signal. Second, while 
we have assumed a specific neutron star mass of 1.4M , Fig.[2T]tells us that /in decreases 
with increasing stellar mass, i.e. heavier neutron stars will emit weaker GWs. Here we re- 
call that from the selected pulsars only the mass of PSR J0437-4715 is known (within some 
accuracy lfT25l[T26l ). The masses of PSR J2124-3358 and PSR J1024-0719 are unknown. 
Third, in the present study we assume a very simple model of stellar matter consisting only 
beta equilibrated nucleons and light leptons (elections and muons). On the other hand, 
in the core of neutron stars conditions are such that other more exotic species of particles 
could readily abound. Such novel phases of matter would soften considerably the EOS of 
stellar medium II 1031 leading to ultimately more compact and gravitationally tightly bound 
objects which could withstand larger deformation forces (and torques). Lastly, the exis- 
tence of quark stars, truly exotic self-bound compact objects, is not excluded from further 
considerations and studies. Such stars would be able to resist huge forces (such as those 
resulting from extremely rapid rotation beyond the Kepler, or mass-shedding, frequency) 
and as a result retain their axial symmetric shapes effectively dumping the gravitational ra- 
diation (e.g. Ref. [4]). At the end, we recall that Eq. (1401 implies that the best possible 
candidates for gravitational radiation (from spinning relativistic stars) are rapidly rotating 
pulsars relatively close to Earth (/in ~ <3>22^ 2 / r )- Increasing rotational frequency (and/or 
decreasing distance to detector, r) would alter the results shown in Figs.|2T]and|22]in favor 
of a detectable signal by the current observational facilities (e.g. LIGO). On the other hand, 
for more realistic and quantitative calculations, the neutron star quadrupole moment must 
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Table 8: Properties of the rapidly rotating pulsars considered in this study. The first column 
identifies the pulsar. The remaining columns exhibit the following quantities: rotational 
frequency; first derivative of the rotational frequency; distance to Earth; corresponding 
reference. 

Pulsar u(Hz) v(Hz s _1 ) r(kpc) Reference 

PSRB 1937+21 641.93 -4.33 x 1(T 14 3^60 01113011 
PSR J1748-2446ad 716.35 - 8.70 HBIfTBH 



be calculated numerically exactly by solving the Einstein field equations for rapidly rotat- 
ing neutron stars. (Such calculations have been reported, for instance, by Laarakkers and 
Piosson IH7I .) 



7.3 Gravitational waves from rapidly rotating neutron stars 

For rapidly rotating neutron stars the estimate of the quadrupole deformation in Eq. (|42l) is 
no longer valid. Moreover, the moment of inertia has to be calculated differently as well. 
Often, one uses the spin-down rate to estimate the strain amplitude for fast pulsars. Provided 
that the spin-down rate, u, for a given pulsar is known, ho could be estimated from 1112811 

^»Y /2 . (44) 



Here we should mention that in such calculations one assumes that the only mechanism 
contributing to the pulsar's observed spin-down is gravitational radiation. However, other 
mechanisms could also account for the star's observed decrease in rotational frequency such 
as magnetic dipole radiation, and particle acceleration in the magnetosphere II129II . Despite 
these uncertainties, calculations of gravitational wave strain amplitude through Eq. (1441 are 
still important because, in addition to providing a rather conservative upper limit on the ex- 
pected gravitational radiation, they also serve to estimate another very uncertain but impor- 
tant quantity - the ellipticity e. So far the fastest pulsars reported are the PSR B 1937+21 Q 
and PSR J1748-2446ad [5]. Their known properties are summarized in Table[8] 

The RNS code ll46ll was used to calculate the principal moment of inertia (and other 
properties) of rapidly rotating neutron stars. The neutron star moment of inertia for the two 
fastest pulsars is shown in Fig|T4l Panel (a) shows the moment of inertia of PSR B 1937+2 1 
and panel (b) displays the moment of inertia of PSR J1748-2446ad. As already observed 
previously, the moment of inertia increases with rotational frequency, while the range of 
possible neutron star configurations decreases (see, for instance, Refs. ||59l 601). 

A particularly interesting example is the neutron star rotating at 642 Hz (3J- Since its 
first observation in 1982, this pulsar has been studied extensively and an observed spin- 
down rate has been measured (see Table 1). Using the spin-down rate, an upper limit on 
the gravitational wave strain amplitude was obtained. The spin-down rate corresponds to a 
loss in kinetic energy at a rate of E = A.ir 2 I zz v\i>\ ~[0.6 - 3.1] x 10 36 erg/ s, depending on 
the EOS. Assuming that the energy loss is completely due to the gravitational radiation, the 
gravitational wave strain amplitude can be calculated through Eq. (l44l . Similar calculations 
for this pulsar and others with an observed spin-down rate have been performed in the 
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Figure 23: (Color online) Gravitational wave strain amplitude tiff (left panel) and ellipticity 
(right panel), deduced from the spin-down rate of PSR B 1937+21. See text for details. 

past 111081 11281 . These calculations have provided estimates for the gravitational wave 
strain amplitude of selected pulsars for which the spin-down rates are known, and also upper 
bounds for their ellipticities using the quadrupole model 111081 . However, such calculations 
simply used the "fiducial" value of 10 45 g cm 2 for the moment of inertia I zz in all estimates. 
On the other hand, the neutron star moment of inertia is sensitive to the details of the 
EOS of stellar matter, and especially to the density dependence of the nuclear symmetry 
energy [60]. Moreover, I zz increases with increasing rotational frequency (see, e.g. Fig.lT5T> 
and the differences with the static values of the moment of inertia could be significant, 
particularly for rapidly rotating neutron stars. 

The gravitational wave strain amplitude of PSR B 1937+21 is shown in Fig. [23] (left 
panel). Because the MDI EOS is constrained by available nuclear laboratory data, our re- 
sults with the x = and x = — 1 EOSs allowed us to place a rather conservative upper limit 
on the gravitational waves to be expected from this pulsar, provided the only mechanism 
accounting for its spin-down rate is gravitational radiation. Under these circumstances, the 
upper limit of the strain amplitude, hff, for neutron star models of 1.4M is in the range 
h s d = [2.24-2.61] xl0~ 27 . Similarly, we have constrained the upper limit of the ellipticity 
of PSR B 1937+21 to be in the range e sd = [2.68 - 3.12] x 10" 9 (Fig.[23 right panel). 

8 Summary and outlook 

In this chapter we have summarized our recent studies on properties of (rapidly) rotating 
neutron stars and the gravitational waves expected from deformed pulsars using several 
nuclear EOSs constrained partially by terrestrial laboratory experiments. In particular, the 
latest heavy-ion reaction experiments have constrained partially the density dependence of 
the nuclear symmetry energy and thus the EOS of neutron-rich nuclear matter. These lim- 
its, while being incomplete and still suffer from some remaining uncertainties, can already 
provide some useful information about the possible stable configurations of (rapidly) rotat- 
ing neutron stars and the gravitational radiation expected from them. Specifically, we have 
studied properties of (rapidly) rotating neutron stars employing four nucleonic EOSs and 
find that the rapid rotation affects the neutron star structure significantly. It increases the 
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maximum possible mass up to ~ 17% and increases/decreases the equatorial/polar radius 
by several kilometers. The neutron star moment of inertia has been studied for both slowly 
and rapidly rotating models within a well established formalism. We found that the moment 
of inertia of PSR J0737-3039A is limited in the range of I = (1.30- 1.63) x 10 45 (# cm 2 ). 
The fractional momenta of inertia AI/I of the neutron star crust are also constrained. It is 
also found that the moment of inertia increases with rotational frequency at a rate strongly 
depending upon the EOS used. Additionally, rotation reduces central density and proton 
fraction in the neutron star core, and depending on the exact stellar mass and rotational 
frequency could effectively close the fast cooling channel in millisecond pulsars. This cir- 
cumstance may have important consequences for both the interpretation of cooling data and 
the thermal evolution modelling. 

We have reported predictions on the upper limit of the strain amplitude of the gravi- 
tational waves to be expected from elliptically deformed pulsars at frequencies < 300Hz. 
Our results are intended to provide guidance to the ground-based gravitational wave ob- 
servatories. By applying an EOS with symmetry energy constrained by recent nuclear 
laboratory data, we obtained an upper limit on the gravitational-wave signal to be ex- 
pected from several pulsars. Depending on details of the EOS, for several millisecond 
pulsars 0.18/cpc to 0.35Kpc from Earth, the maximal ho is found to be in the range of 
~ [0.4 - 1.5] x 10~ 24 . Finally, from the spin-down rate of PSR B1937+21 we have de- 
duced the upper limit of the strain amplitude, hff, for neutron star models of 1.4M to be 
in the range hff = [2.24 — 2.61] x 10~ 27 . We have also constrained the upper limit of the 
ellipticity of PSR B 1937+21 to be in the range e sd = [2.68 - 3.12] x 10~ 9 . These pre- 
dictions set the first direct nuclear constraints on the gravitational waves from elliptically 
deformed pulsars. 

Looking forward, new experiments in terrestrial nuclear laboratories are expected to 
improve our understanding about the EOS of dense neutron-rich nuclear matter dramatically 
in the next few years. The EOS of neutron-rich matter will be better constrained in a wide 
density range. Most of the studies presented in this review will then be further refined. 
Conversely, rapid progress in astrophysical observations of neutron stars and gravitational 
waves surely will also help us better understand the EOS of neutron-rich matter, and thus 
stimulate the progress in nuclear physics. Eventually, the ultimate goal of understanding 
thoroughly all mysteries of pulsars can only be realized by utilizing the progress made in 
both nuclear physics and astrophysics. 
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